Theromelectricity in Graphene: Effects of a gap and magnetic fields 
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We calculate the thermopower of monolayer graphene in various circumstances. First we show that 
experiments on the thermopower of graphene can be understood quantitatively with a very simple 
model of screening in the semiclassical limit. We can calculate the energy dependent scattering 
time for this model exactly. We then consider acoustic phonon scattering which might be the 
operative scattering mechanism in free standing films, and predict that the thermopower will be 
linear in any induced gap in the system. Further, the thermopower peaks at the same value of 
chemical potential (tunable by gate voltage) independent of the gap. Finally, we show that in the 
semiclassical approximation, the thermopower in a magnetic field saturates at high field to a value 
which can be calculated exactly and is independent of the details of the scattering. This effect might 
be observable experimentally. 

PACS numbers: 72.80.Vp, 79.10.N-, 72.15Jf 
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Graphene has become one of the most studied elec- 
tronic systems in recent years. The interest in this sys- 
tem stems from the existence of a relativistic band struc- 
ture which gives rise to several interesting electrical prop- 
erties Very high values of electron mobility com- 
bined with robust mechanical properties make graphene 
a promising candidate for use in several devices and tech- 
nologies. Two recent developments have proven very sig- 
nificant in this regard: The ability to introduce a tunable 
gap in single sheets of graphene @, H| and the fabrica- 
tion of free standing sheets The former, in principle 
enables one to fabricate electronic devices using single 
sheets of graphene while the latter is an important step 
towards increasing the mobility of electrons by eliminat- 
ing scattering effects from the substrate. 

An enhanced mobility (and hence electrical conductiv- 
ity) is important for thermoelectric applications as well. 
In addition to the technological significance of the phe- 
nomenon of thermoelectricity, the thermopower of a sys- 
tem is often used to shed light on various properties of 
the elementary charge carriers. These include the sign 
of their charge, the principle scattering mechanism and 
in certain limits, their entropy. Thus, the thermopower 
in graphene has also been studied extensively in the last 
few years both experimentally [H-0] and theoretically 0- 
131 ] . Interesting behavior such as the violations of the 
Mott formula at high temperatures have been found and 
attributed to the non-degenerate behavior of carriers. 
Most theoretical calculations have focussed on monolayer 
graphene without a gap with screened impurities being 
the dominant source of scattering. This is indeed the sit- 
uation for graphene on substrates of SiC>2 on which most 
experiments have been carried out. However, with the 
two new developments mentioned in the previous para- 
graph, it has become important to consider the effect of 
a gap and other scattering mechanisms, which is one of 
the objectives of this paper. 

A magnetic field can have interesting effects on thermo- 



electricity. It can produce oscillations @, EH or steps 14 1 
in the thermopower under suitable circumstances which 
can provide information about the energy levels of the 
carriers of the system. The effect of magnetic fields on the 
thermopower of graphene has also been studied, mainly 
at low carrier densities, where the cyclotron frequency 
is high and oscillations have been seen indicative of the 
formation of Landau levels 0-0, E3- However, not much 
attention has been paid to the effect of a magnetic field 
at higher values of filling when the cyclotron frequency 
is low and semiclassical physics can apply at reasonable 
temperatures. We show that in such a situation the ther- 
mopower can saturate at high field, an effect that might 
be observable experimentally 

In this paper, we consider semiclassical transport in 
graphene taking into account the effects of a gap, various 
scattering mechanisms and magnetic fields. We first show 
that it is possible to obtain quantitative agreement with 
thermopower measurements on graphene with a simple 
model of scattering in the semiclassical limit. While cal- 
culations along the same lines have been performed ear- 
lier our calculation has the virtue of simplicity 
which allows us to calculate the scattering rate exactly 
and consequently make a fit to experimental data using 
only a weakly temperature dependent screening length 
and no other parameter. We then move on to consider the 
case of scattering due to acoustic phonons, which might 
be the dominant scattering mechanism in free standing 
films due to the absence of impurity effects from the sub- 
strate. Here, we find predict that the thermopower is 
proportional to the gap and obtain an exact expression 
in the limit of small gaps. This is quite distinct from 
what happens in the presence of substrates and might be 
used to characterize the gap. Finally, we calculate the 
effect of magnetic fields and find that in the semiclassi- 
cal limit, the thermopower saturates to a value at high 
field, that can be calculated exactly, independent of the 
details of scattering. We show that for reasonable values 
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of various parameters, this saturation can occur at a few 
Tesla of field, making it relevant experimentally. 

Graphene has a hexagonal Brillouin zone. At the six 
corners of the zone, the dispersion is given by E(k) = 
±?ii>/|k| where Vf is the fermi velocity. The density of 
states is 
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In linear response response theory in the semiclassical ap- 
proximation, the transport coefficients relevant to ther- 
moelectric response at temperature T are given by: 
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where Vi is the i component of the velocity, p, the 
chemical potential and / the Fermi function. The ther- 
mopower 
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The momentum integrals are carried out in the vicinity of 
all six corners of the Brillouin zone. The carrier density 
p can be related to p through the expression 
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-f(E)]g(E)dE. (4) 



p is the excess carrier concentration over the completely 
filled lower band and p = when p = 0. 

As mentioned earlier, for graphene films on SiC>2 sub- 
strates the dominant scattering mechanism is due to 
charged impurities screened by the conduction electrons. 
While it is possible to obtain the screened dielectric func- 
tion using techniques like the random phase approxima- 
tion [8l4l3j , we will assume a simple screened potential of 
the form 
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with a screening length £. As will be seen, a weakly 
temperature dependent £ enables us to make quantita- 
tive fits to the experimental data. The main advantage 
of this potential is that the relaxation time can be calcu- 
lated exactly using the Fermi golden rule yielding a fairly 
simple form 
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where q is the charge of the impurity, pi is the impurity 
concentration and E(k) is given by Eqn.[SJ When £ = oo, 
one obtains the known result t(E) oc \E\ (TBI ]. 
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FIG. 1: Thermopower of graphene as a function of chemi- 
cal potential at different temperatures due to the scattering 
potential of Eqn. [5] The calculated peak values of the ther- 
mopower (diamonds) at every temperature can be directly 
compared to the experimental data (asterisks) of Ref. 0] (in- 
set). The only fitting parameter is the screening length £. 
The values £ = 8.5 nm for T = 80 - 160 K and £ = lOnm for 
T — 180 — 300 K have been used here to obtain best fits to 
the experimental data. 



Eqn. [6] can be used to calculate the thermopower. The 
required integrals have to be evaluated numerically and 
the thermopower as a function of p at different temper- 
atures is show in Fig. [TJ The only parameter that is re- 
quired to obtain this data is the screening length £, which 
we consider a fitting parameter. The thermopower has 
experimentally been measured as a function of gate volt- 
age 0] , which can be converted to p only from knowledge 
of the sample dependent capacitance. However, the peak 
thermopower at every temperature can be read off di- 
rectly and compared to experiments. We find very good 
agreement over a fairly large range of temperatures using 
a screening length of 8.5 nm for 80-160 K and 10 nm for 
180-300 K. 

The thermopower due to this screened potential can 
also be calculated in the presence of a gap A and does 
not show a very strong dependence on it at least for 
A <C ksT. A much more interesting situation arises 
when the dominant scattering mechanism is acoustic 
phonon scattering. This mechanism is likely to be rel- 
evant for free-standing graphene films of the sort that 
have been fabricated recently. It can be shown from the 
Fermi golden rule that the relaxation time t(E) cx l/\E\ 
for this mechanism. Incidentally, this is also the energy 
dependence of t(E) for localized impurity scattering. In 
the absence of a gap, the energy dependence of the prod- 
uct of the density of states with the square of the velocity 
in the integrals for is exactly cancelled by that com- 
ing t(E). Thus, 
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since the integrand is an odd function of E — fi. The ther- 
mopower is identically zero with acoustic phonon scat- 
tering independent of the chemical potential (and hence 
carrier density) fisj] . To demonstrate the effect of a gap, 
we assume dispersion of the form 



E(k) = ±Jh v 2 \k\ 2 + A 2 . 
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The magnitude of the velocity now picks up an energy 
dependence 
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g(E) is still given by Eqn. Q] and p can be calculated in 
terms of fi to yield 
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where Li2(x) is the second polylogarithmic function. 

The integrand in the expression for L X2 X is no longer an 
odd function of E — (j,. Further, the limits of integration 
also change. However, the form of t{E) stays the same 
for acoustic phonon scattering due to g(E) not changing. 
The resultant integration gives a non-zero value of the 
thermopower. This value can be obtained analytically in 
the limit of small gaps, A <C /i and A <C k B T and is 
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roughly linear to reasonably large values of the gap. The 
thermopower as a function of chemical potential is shown 
in Fig. |3l A curious fact is that the maximum appears 
at a value of \ijk B T = 1.54 roughly independent of the 
band gap. This can be seen directly for very low band 
gaps from Eqn. [TT] but holds more generally for larger 
gaps as well. The maximum value of the thermopower 
can be seen to be of the order of k B /e for sufficiently large 
values of the gap and is shown in the inset of Fig. [3] It 
can also be seen that this peak thermopower is also in- 
ear in the gap up to far large values of the gap. Thus, 
the thermoelectric properties of a gapped free standing 
graphenc sheet arc likely to be quite different from those 
of a sheet on a substrate. In particular, the linear rela- 
tion between the gap and the peak thermopower in the 
former case might provide a method to estimate the size 
of the gap. 
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FIG. 2: Thermopower as a function of band gap for acoustic 
phonon scattering. The different curves correspond to differ- 
ent values of chemical potential. 

The thermopower can be calculated numerically for 
larger values of the gap. The results are shown in 
Fig- El from which it can be seen that the thermopower is 



FIG. 3: Thermopower for different values of chemical poten- 
tial for acoustic phonon scattering. The different curves cor- 
respond to different values of the gap. It can be seen that the 
peak occurs at approximately /i = 1.5fcsT independent of the 
value of the gap. (Inset) The peak value of the thermopower 
as a function of the gap. 

We now focus our attention to the case of magnetothcr- 
mal transport. The thermopower of graphene has been 
measured in the presence of a magnetic field at low tem- 
perature where oscillations have been observed as a func- 
tion of the field consistent with quantization of energy 
levels @, 0- Here we focus on semiclassical transport 
in the presence of a magnetic field which is relevant to 
experiments at higher temperature. The field B is as- 
sumed to be perpendicular to the plane of the sheet. We 
consider the effect of magnetic fields in the absence of a 
gap. The solution of the semiclassical equations of mo- 
tion gives orbits in which the energy of the electron is a 
constant. The linear response formulae Eqn. [5]then have 
to be modified to include the average of the velocities 
along these orbits. Operationally, this involves replacing 
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one of the velocit y o p erators by its average value in an 
orbit given by [lg, [l7[ 
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where t{E) is the relaxation time. It can be shown that 
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FIG. 4: Thermopower as a function of magnetic field for 
T = 273 K and h/UbT = 1.59, for acoustic phonon scat- 
tering and screened scattering (Eqn.[5|. For acoustic phonon 
scattering, the deformation potential for graphene was taken 
to be 18 eV, the velocity of LA phonons 2.1 x 10 6 cm/s and 
mass density 7.6 x 10 -7 Kg/m 2 . For screened impurity scat- 
tering, the impurities were assumed to have a charge of 2e 
and concentration of 1 ppm, and screening length of 10 nm. 

the following relations result for the various transport 
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coefficients when t(E) explicitly depends on |-E7| 
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where j3 = ewji?[20|. Note that while the integrals over 
the two bands for the diagonal coefficients add, the 
ones for the off-diagonal coefficients subtract. This 
is due to the fact that the sense of rotation of the or- 
bits is opposite in the two bands. We have verified that 
the above integrals for screened impurity scattering yield 
values of conductivity and hall resistance as functions of 
fi that show very good qualitative agreement with the 
experimental data of Ref. The Onsagcr relations 
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The thermopower in the presence of a magnetic field is 
given by the formula 
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An analytic expression can be obtained for the saturation 
value of the thermopower at large magnetic fields. 
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where 1^3(2) is the third poly logarithmic function. For 
different scattering mechanisms, the field dependence of 
the thermopower has to be calculated numerically. This 
is shown for two scattering mechanisms, acoustic phonon 
scattering and screened impurity scattering in Fig. [4] for 
realistic values of the relevant parameters. The ther- 
mopower approaches the value given by Eqn.[l5]for large 
enough fields in both cases. The saturation value is inde- 
pendent of the scattering mechanism as can be seen from 
Eqns.[T31 For sufficiently large fields, (3/E dominates the 
factor 1/t(E) where the value of the integrand is appre- 
ciable. In this limit, S ~ Lj. 2 /L^, which can be shown 
to be equal to Eqn. [151 Effects of the Zeeman splitting 
of the electrons can be neglected even up to fields of 10 
T since the corresponding energy is still much smaller 
compared to the Fermi energy. For the parameters used, 



the saturation field is of the order of a few Tesla as can 
be seen from Fig. SI which could make the effect experi- 
mentally observable. 
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